Using Full Information When Computing Modes of Post-Newtonian Waveforms Prom 
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The increasing sophistication and accuracy of numerical simulations of compact binaries (espe- 
cially binary black holes) presents the opportunity to test the regime in which post-Newtonian (PN) 
predictions for the emitted gravitational waves are accurate. In order to confront numerical re- 
sults with those of post-Newtonian theory, it is convenient to compare multipolar decompositions 
of the two waveforms. It is pointed out here that the individual modes can be computed to higher 
post-Newtonian order by examining the radiative multipole moments of the system, rather than by 
decomposing the 2.5PN polarization waveforms. In particular, the dominant (I = 2,m = ±2) mode 
can be computed to 3PN order. Individual modes are computed to as high a post-Newtonian order 
as possible given previous post-Newtonian results. 

PACS numbers: 04.25.Nx, 04.25.Dm 



I. INTRODUCTION 

The first generation of laser interferometric gravita- 
tional wave detectors is now operating at or near their 
design sensitivities [H, [3, H, 0, Q ■ One of the most promis- 
ing sources that they may detect are the inspiral and 
merger of compact binaries systems with black holes or 
neutron stars. One of the primary means of detecting and 
interpreting the signals from such systems will be the use 
of a matched filtering technique, which requires accurate 
templates with which to match a theoretical model to a 
noisy signal. Until the last several orbits prior to merger, 
the post-Newtonian (PN) approximation is expected to 
be sufficient to generate accurate templates, while numer- 
ical relativity will be needed to help construct accurate 
templates covering the late inspiral and merger phases. 

Estimates of the accuracy of post-Newtonian templates 
and the effect of this accuracy on detecting gravita- 
tional waves have led the post-Newtonian expansion to 
be pushed to high order 0, 0, H, • Currently the equa- 
tions of motion and the gravitational wave energy flux 
are known through 3.5PN order and the polar- 

ization waveform amplitudes are known through 2.5PN 
order [l^, As the post-Newtonian expansion is a 

slow-motion, weak-field approximation to general rela- 
tivity, and the motion is not so slow and the fields not so 
weak during the last orbits prior to merger, it is legiti- 
mate to ask where the post-Newtonian expansion breaks 
down. 

Recent breakthroughs in numerical relativity [3, fisl . 
[igI . \vf\ have finally provided the possibility of testing the 
validity of the post-Newtonian expansion in the late in- 
spiral. Initial studies focused on qualitative comparisons 
between post-Newtonian generated waveforms and nu- 
merical simulations of the final two to four orbits prior 
to merger El) ■ More recently, attempts have been 
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made to quantify the difference between post-Newtonian 
generated waveforms and numerical simulations of a non- 
spinning equal-mass binary inspiral lasting more than 
eight orbits [l!], [l^, ■ As long inspiral simulations are 
computationally expensive to perform, it is important 
to quantify the errors in the post-Newtonian generated 
waveforms and determine how they influence the detec- 
tion of gravitational waves. One would also like to know 
the accuracy with which the parameters such as the indi- 
vidual masses and spins of the binary can be determined 
from the observed waveform. 

One of the principal ways of quantifying the accuracy 
of the post-Newtonian approximation is to compare the 
amplitude and phase of the gravitational waveform with 
that from numerical simulations. Unlike gravitational 
wave detection which is sensitive only to the waveform in 
one particular direction, comparisons between numerical 
simulations and post-Newtonian waveforms can be per- 
formed over the entire sphere. Therefore it is convenient 
to project the waveforms onto spin-weighted spherical 
harmonics and compare the individual components. Pre- 
vious studies [H, [ill, however, have computed the spin- 
weighted spherical harmonic components of the post- 
Newtonian amplitude in such a way that has not made 
full use of current post-Newtonian results. 

This situation has arisen primarily because of how 
the post-Newtonian gravitational wave amplitude is pre- 
sented. Often the final result of the post-Newtonian 
computation of the waveform is given as the polariza- 
tion waveforms as a function of the direction to the ob- 
server (as this is most useful when constructing tem- 
plates) . Currently the polarization waveforms are known 
to 2.5PN order in amplitude [12, [III- previous com- 
parisons between post-Newtonian and numerical wave- 
forms, it was the polarization waveforms that were pro- 
jected into spin-weighted spherical harmonics in order 
to determine the post-Newtonian components. The 
post-Newtonian polarization waveforms are computed by 
truncating an infinite sum of radiative multipole mo- 
ments at a given post-Newtonian order [23|. As will 
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be demonstrated in this paper, it makes more sense to 
project the waveform onto spin-weighted spherical har- 
monics prior to truncating the series as this retains the 
fuU information currently known. In particular, the dom- 
inant {I — 2, m — 2) component of the waveform can be 
computed to 3PN order. 

While knowledge of a particular spin-weighted compo- 
nents of a gravitational waveform is of limited value for 
computing templates for detection and characterization 
of gravitational waves, it can be used to quantify the 
difference between post-Newtonian and numerical wave- 
forms. In addition, if the 3PN correction to the am- 
plitude of the {I — 2,m = 2) component significantly 
improves the agreement between the post-Newtonian 
and numerical waveforms, this would suggest it might 
be worth the effort needed to compute the polarization 
waveforms to 3PN order. It has been shown that higher 
post-Newtonian corrections in the amplitude can improve 
detection rates i25i. i26i. i27i] as well as improve parameter 
extraction [ll, [illSjSlUll [sl, 113 . 

The relationship between the gravitational waveform 
and radiative multipole moments is reviewed in Sec. [HI 
It is pointed out that the the spin-weighted spherical har- 
monic components of the waveform can be computed di- 
rectly from the radiative mulitpole moments of the sys- 
tem in Sec. Ill B[ Current post-Newtonian results that af- 
fect the computation of the radiative multipole moments 
are summarized in Sec. IIIII The radiative multipole mo- 
ments are computed to as high a post-Newtonian order as 
possible in Sec.|TVl Conclusions are presented in Sec.lVl 



II. GRAVITATIONAL WAVEFORM AND 
SPIN- WEIGHTED SPHERICAL HARMONICS 

A. Notation 

Let = {cT, X, Y, Z) be coordinates in an asymptot- 
ically flat radiative coordinate system, with eV the unit 
timelike normal, and {ex , eV, ez) the spatial orthonormal 
coordinate basis vectors. Let (cT, R, 8, $) be the corre- 
sponding spherical coordinate system with corresponding 
basis vectors (br, ee, e$). Let Tfj — T ~ R/c denote re- 
tarded time in radiative coordinates. 

Greek letters are used for spacetime indices, Latin let- 
ters for spatial indices. As indices can be raised or low- 
ered using the Minkowski metric, components of tensors 
will typically have down indices. Let Ni be a component 
of the unit radial vector cr. A capital letter for an index 
denotes a multi-index (i.e. Tl = Ti^i^...i^). For a vector 
V , let Vl be a product of components of the vector (i.e. 
Vl = Vi^Vi^ ■ ■ ■ Vig) . Repeated spatial indices (includ- 
ing multi-indices) are to be summed over. Tensors that 
are fully symmetric and trace-free (STF) will be denoted 
with capital script letters (e.g. Xy). Symmetrization, 
anti-symmetrization, and STF projection of indices are 
denoted by (), [], and <> respectively, with underlined 
indices between the delimiters to be excluded from the 



operation (e.g. 2T(^abc) = Tabc + Tcba)- 

Let (fc, /, rh) be a set of null vectors defined by: 

k = -^{cT + en), (1) 

I = 7f(e^T-ei?), (2) 

rh ~ — (ee + ie$) . (3) 

The spin-weighted spherical harmonics are defined in 
terms of the Wigner d-functions by 



,y'^"(e,$) = (-i)^^/^d^^(e)e™*, (4) 



where 



dL(0) =^{l + m)\{i-m)\{l + s)\{l-s)\ 
^ ^ (-l)''-(sinf )2'''+"-™(cosf )2^+™-^-2fe 



k=ki 



k\{i + m- ky.{i - s - fc)!(s - m + fc)! 



,(5) 



where fc; = max(0, m — s) and kf = min(^ -I- m, £ — s). 

B. Polarization Waveforms and Modes 

The asymptotic waveform hfj^ can be decomposed into 
two sets of symmetric trace-free (STF) radiative multi- 
pole moments as [23 | 

AG ^ r 1 

i=2 

2£ 1 

Here IAl{Tr) are the mass- type moments and Vl{Tr) 
are the current-type moments. In Sec. IIII Al these ra- 
diative multipole moments will be related to multipole 
moments describing the source in the near-zone. Note 
that higher multipole moments contribute to the wave- 
form at higher post-Newtonian order, so that at any finite 
post-Newtonian order, only a finite number of multipoles 
contribute. The transverse-traceless (TT) projection op- 
erator liijmn is given by 



IT — P P 



where Pij — 6ij — NiNj. 



1 



- P P 



(7) 



Given an orthonormal triad (A^, P, Q), the polarization 
waveforms can be given by 



TT 
mn ■ 



(8) 
(9) 
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A natural (but by no means unique^) choice for the triad 
is P = ee and Q = e$. It is then straightforward to 
show that 



(10) 



where * denotes complex conjugation. It will now be 
shown how h^ — ihx can be decomposed into modes using 
spin- weighted spherical harmonics of weight -2 

cx; t 

h+-ihx^Y, /i'"-2>"'™ (e, $) . (11) 

An alternative expression for the waveform is given 
by HI 



~ c-^R 



f — O ^y^ — — ^ 



E2,lm 



Im irp_\rriB2,lm 



(12) 



where and L^^ are pure-spm tensor harmon- 

ics, and where the mass multipole moments [/^"'(T/j) and 
current multipole moments V^™'{Tfi) are related to their 
STF counterparts by [H] 



16. ^^3^ 



(2£-hl)!!y 2^(^-1) 
-327rf / 



{1 + 2) 



(2£-H)!!y 2^(£-|-l)(£- 1) 



v^yr*. (14) 



where 3^^™ are the STF spherical harmonics which are 
related to the scalar spherical harmonics by 

r^"'(e,$) ^j^i^A^L- (15) 

The pure-spin tensor harmonics are related to the spin- 
weighted spherical harmonics by (2^ 1 



1 



B21m ~i 



(_2r^™™,m, + 2Y'"'m*m*) , (16) 



^ {-2Y'"^m,m, ~ ^r'-m*™;) . (17) 



Combining Eqs. dTO]), (HH), and ^ yields 

-2 m——t ^ 



:^i-^'"(Tfl)|_2r^"(e,<i>)(i8) 



Thus the spin-weighted spherical harmonic components 
of the waveform are given by 



;^^(t/^™(T.)-V-(T,)). (19) 



Note that the spin-weighted spherical harmonic com- 
ponents can be computed directly from the radiative mul- 
tipole moments. It is not necessary to compute the wave- 
form as a function of position, and then project out the 
components as is commonly done in the recent literature. 
In fact, not only are these extra steps unnecessary, they 
cause individual components to be truncated at a lower 
post-Newtonian order than they need to be. In particu- 
lar, the dominant h^^ component is truncated to 2.5PN 
order, when it can be computed (as will be shown below) 
to 3PN order directly from the radiative mass quadrupole 
moment. 

This situation arises as the complete polarization wave- 
forms /14. and /ix are known only through 2.5PN or- 
der. But as discussed in [l^], the computation of the 
3PN waveform is not limited by the post-Newtonian or- 
der of the dominant mass quadrupole contribution, but 
by the order of higher multipoles. Therefore to obtain 
the spin- weighted spherical harmonic components of the 
waveform to as high a post-Newtonian order as possi- 
ble, they should be computed directly from the radia- 
tive multipole moments. As discussed in Sec. |V] below, 
this will allow post-Newtonian results to be probed to a 
higher order when comparing with numerical simulations. 
Table U summarizes the post-Newtonian order to which 
each radiative multipole moment needs to computed to 
obtain the full 2.5PN order polarization waveforms, and 
the post-Newtonian order to which they are currently 
known. 



III. SUMMARY OF CURRENT 
POST-NEWTONIAN RESULTS 

The post-Newtonian approximation is a slow-motion, 
weak-field approximation to general relativity with an 
expansion parameter e ~ (v/c)'^ ^ (Gm/rc^). For a 
binary system of two point masses mi and m2, v is 
the magnitude of the relative velocity, m = mi + 1712, 
and r is the separation. In order to produce a post- 
Newtonian waveform, it is necessary to solve both the 
post-Newtonian equations of motion describing the bi- 
nary, and the post-Newtonian equations describing the 
generation of gravitational waves. For a complete review 
of post-Newtonian methods applied to inspiralling com- 
pact binaries, see [35j . 



Post-Newtonian Computation of Radiative 
Multipoles 



^ The 2.5PN polarization waveforms in [T^ are defined with P = 
—e.^ and Q = eg evaluated at © = j, <& = 7r/2. 



In Sec. Ill BI the spin-weighted spherical harmonic com- 
ponents have been related to the radiative multipole mo- 
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TABLE I: The post-Newtonian order to which each radiative 
multipole moment must be computed to obtain the full 2.5PN 
polarization waveforms h+ and h^, and the order to which 
they are known. 

Radiative Multipole PN order for /i+,x Known FN order 
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3 


*Jij 
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1'ijk 
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<!Jijk 


1.5 


1.5 


"I-ijkl 


1.5 
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J'ijkl 
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"I-ijklm 


1 


1 




0.5 


1 


■^ijklmn 
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<!Jijklmn 





1 


-^ijklmno 





1 


Il{1>7) 




1 


jL{i > 6) 




1 



"The radiative current octupole could be computed to 2PN or- 
der by a computation of the 2PN correction to the source current 
octupole. 

ments. In order to use Eq. it is necessary to relate 
the radiative moments to a description of the compact bi- 
nary system. The post- Newtonian wave generation for- 
malism has been developed [sl, [sS, S iHi, SO, Sll, SI 
in a systematic manner to relate the radiative multipole 



moments {Wl,Vl} to a set of six STF source moments 
{Tl, JL,y^L,XL,yL, Zl}, which can be computed from 
the stress-energy pseudo-tensor of the material and grav- 
itational fields of the source. A set of two canonical 
source moments {Ml,Sl} can be computed as a non- 
linear functional of the source moments as an intermedi- 
ate step between the source moments and the radiative 
moments. The radiative moments are then given as non- 
linear (and even nonlocal) functionals of the canonical 
moments. It turns out that two of the source moments, 
the source mass moments Ii, and the source current mo- 
ments J7l are dominant, while the other four parameter- 
ize a gauge transformation and only make a contribution 
to the canonical source moments starting at 2.5PN or- 
der. To the post-Newtonian order needed in this paper, 
only the canonical mass quadrupole moment contains a 
correction from its corresponding source moment [43l . 

M., = J., + ^ (W^'h,, - W(i)xg^) + Oie'/'), (20) 

where Tjf^ denotes the pth time derivative of 7^ and 
0(6*^) denotes the fcth-ordcr and higher post-Newtonian 
corrections that are unknown. The other canonical 
source moments are simply related to the source moments 

Ml = lL + 0(e'/'), (21) 
Sl = JL + 0(e5/2). (22) 
The radiative mass quadrupole is given by 
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124627 
44100 



+ 0(e^/2). 



(23) 



The first integral in the above expression is the domi- 
nant radiation tail at 1.5PN order obtained in [s^. The 
2.5PN non-linear memory integral has been obtained 
in [sl, [13, m, H^, 113. The other non-linear interac- 
tions at 2.5PN order were obtained in |4(]||. F inallv the 
3PN "tail of tail" integral was derived in |4l| . The tail 
integrals involve nonlinear interactions with the mass 
monopole M (equivalent to the ADM mass) of the sys- 
tem. The tail integrals also contain a freely-specifiable 
constant tq that corresponds to the choice of the origin 
of radiative time T with respect to harmonic time t, and 



enters the relation between the retarded time in radiative 
coordinates and the retarded time t — r/c in harmonic 
coordinates (the coordinates in which the equations of 
motion are given) 0, \^ : 



r 2GM T 

TR = t ^ In — 

c \ ro 



(24) 



The remaining radiative multipole moments are given 
by [Si, Sim 
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+ 3 



(25) 

+ 0(e5/2), (26) 

(27) 

+ 0(6^/2), (28) 



3Sl ■ 



fc)b 



25/,.M 



(4) 



0(e 



5/2N 



(29) 
(30) 



r 



B. Adiabatic inspiral of quasi-circular orbits 

Solving the equations of motion yields explicit expres- 
sions for the accelerations of each body in terms of the 
positions and velocities of the two bodies [13, HH, [13, [H, 
0, S S 113, S [si, 113]. The two-body equations of 
motion can then be reduced to relative equations of mo- 
tion in the center-of-mass frame in terms of the relative 
position X and velocity v [6l|. The relative acceleration 
a is currently known through 3.5PN order [13, [llj. The 
effects of radiation reaction (due to the emission of grav- 
itational waves) enters the relative acceleration starting 
at 2.5PN order. 

The emission of gravitational radiation causes the or- 
bits of an isolated binary system to circularize [621] . Thus 
for the remainder of this paper the orbital evolution of 
the binary will be modeled as a slow adiabatic inspiral 
of a quasi-circular orbit. In addition, it will be assumed 
that the individual compact objects are non-spinning. 

The orbital plane is chosen to be the X-Y plane with 
the orbital phase defining the direction of the unit vec- 
tor n = xjr along the relative separation vector by 



cos ex + sm tp ey ■ 



(31) 



Then the motion of the binary can be described using the 
rotating orthonormal triad (n, X,ez) with \ = ez y~ n. 

The relative position, velocity, and acceleration are 
given by 



X = rn, 

V = fn + rujX, 

,2 



a = (f — ruj )n + (rcO + 2fcLi)A, 



(32) 
(33) 
(34) 



where the orbital frequency lu — (f>. Through 2PN order, 
it is possible to model the motion of the binary as a 



circular orbit with the solution r = r ~ uj ~ and rcu'^ = 
—n ■ a. 

At 2.5PN order, however, the inspiral motion must be 
taken into account. The leading order contribution to 
the inspiral of a quasi-circular orbit can be obtained by 
examining the Newtonian orbital energy of a circular or- 
bit, 



E = — V- 
2 



0{e) 



(35) 



where v = mim2/ni? , and the leading order gravitational 
luminosity from a circular orbit. 
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+ 0(e). 



(36) 



and assuming that the energy radiated by the gravita- 
tional waves is balanced by the change in the orbital en- 
ergy (i.e. dE/dt = -L). Then 



. ,dE ,dE, 



64 G^m^ 7/2, 
5 r^d' 



dt dr 

and similarly the orbital frequency changes by 

1/2 



(37) 



^duo ^dr ^ 96 
dr dt 5 



+ 0(e^/2). (38) 



Substituting ^ and ^ into Eqs. ^ and ^ and 
noting that f = O(e^), the following expressions for the 
3PN inspiral relative velocity and relative acceleration 
are obtained: 



64 G^ 



V — rujX h' 



-LU^X - 



32 G^m^ 



n + 0(e^/2), (39) 
0(e^/2). (40) 
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where the 3PN orbital angular frequency is [5J, |5E 



Inverting pTj) order by order, yields 



Gm f , , o / 41 o 



+ T 



-10 



75707 41 , 

\ TT^ + 22 In 

840 64 



(41) 



x<l + x[l u] + [1 



2203 41 2 22 / r 
2520 " 192 - — n 



12 
22 

y 



(44) 



where the post-Newtonian parameter 
7^ 



Gm 



(42) 



and Tq is a gauge constant. 

As 7 is defined with respect to the harmonic-coordinate 
separation r, it is convenient to introduce the frequency- 
related post-Newtonian parameter 



/ Gmuj\ 



2/3 



(43) 



C. Source Multipole Moments 

Defining 6m — mi — TO2, the source mass multipole 
moments for circular orbits are given by [3, [gl, [G^I 



J 

I, 
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0, (46) 
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0(e2), (50) 
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through £ = 8. In Appendix [Al an expression is given for the source mass moments to IPN order for arbitrary £. 
The source current muhipole moments for circular orbits are (sol . 0, [6^ 
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where Li = vmciabXaVb- In Appendix |B] an expression 
is given for the source current moments to leading order 
for arbitrary £. 

Finally, the gauge monopolar moment W for a circular 
orbit is 



1 



(62) 



which is proportional to f ~ 0(e^/^). As W is already a 
2.5PN correction to the source mass quadrupole, it gives 
no contribution in the present case, and the canonical 
moments will simply be given by the source moments. 



A. Evaluation of hereditary terms 



The hereditary terms in the radiative multipole mo- 
ments involve integrals over the entire past history of the 
binary. These integrals fall into two types, the tail inte- 
grals (including the "tail of tail" term) that have loga- 
rithmic terms and physically correspond to the backscat- 
tering of the gravitational waves off the background cur- 
vature 38] ; and the memory integrals that can physically 
be thought of as the re-radiation of the stress-energy of 
the propagating waves [i^, \^ 113] . 

Inserting the canonical multipole moments into the tail 
integrals of the radiative multipole moments yields ex- 
pressions of the form 



IV. RESULTS 

The computation of the spin-weighted spherical har- 
monic components is now an exercise in algebra. The 
evaluation of the tail, memory, and "tail of tail" integrals 
requires special treatment that is described in Sec. IIV Al 
The spin-weighted spherical harmonic components are 
listed in Sec. IIVBI They will contain terms that depend 
upon the freely-specifiable constant rg (see the discus- 
sion between Eqs. and ([211) )• These terms can be 
absorbed into a redefinition of the phase as will be dis- 
cussed in Sec. IIV CI 



dTF{TR -T)e 



-ikLj{Tpt-T) 



ln(^ 

2ro 



(63) 

where k is an integer (the index m of the /i^™ being com- 
puted), 6 is a rational number, and F(T]i — t) represents 
a function whose time dependence enters only through 
its dependence on the orbital frequency. In jH, SI] , it 
has been shown that the oscillatory term in the inte- 
grand combined with the slow adiabatic evolution of the 
orbital frequency allow the integral to be approximated 
by using the orbital frequency of a fixed circular orbit at 
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the current value of Tr, so that 
Ii{Tr) = F(T«)e-'^'=^ 

+ o(ehiO} 



1 y °° dre*^'"" 




+ b 









(64) 



where ^ ^/^'^, is the adiabatic parameter describing 
the decay of the orbit which is 0(e^/^). Equation 
can be evaluated using the identity 



dr In' 



= (7s + In cr) . 



(65) 



where •je is Euler's constant. Using ((65)l to evaluate (|64| 
yields 



/r(r«) ^ -^^(T.)e--^{| 



7B + In f ° I - 6 



(66) 



For reasons that will be explained in Sec. IIV C[ it is con- 
venient to express the above result as 



kuj 



2 \xo 



1 ^ 11 u 



where 



11 



" 2 4, „ 2/Gm 
Inxo = — - -7s - - ln2 + - In ^— 
18 3 3 3 Vc^ro 



,(67) 



(68) 



A similar argument holds for the "tail of tail" integral 
in Eq. ((23)) where a term of the form 



/>oo 

hiTB.) = / dTF{TR-T)e 
Jo 



-ikuj{Tji-T) 



In 



CT 

2^ 



(69) 

is found. In this case the integral can be evaluated with 
the aid of the identity 



(iT(lnT)^e 



to yield 



kuj 



F{TR)e- 



— + (7B+lna)2 
6 



(70) 







r X \ 
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+ ln-H 
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+ ln-H 
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(71) 



The memory integral in Eq. (f26|) leads to an integral 
of the form 



/ F{TR-T)e-''"^('^''-^^dT. (72) 

Jo I— -T 

Using similar arguments as for the tail integrals, Ref. 
has shown that this integral can be approximated by 



kuj 



F{TR)e-'^^ [1 + 0(0] 



(73) 



The memory integrals also lead to an integral of the 
form 



h{TR 



X^{Tr - T)dT. 



(74) 



Unlike the other integrals, it does not depend upon the 
orbital phase, and thus is much more sensitive to the past 
history of the binary. This is the non-linear memory 
effect described in jisl l46l |47| . This integral has been 
evaluated by [H, |46| using the evolution of the frequency 
parameter x found by integrating Eq. 



-1/4 



[l + 0{e)] (75) 



where tc denotes the time of coalescence. With this 
model, 



HTr) 



5 Gm 
64 uc^ 



'[1 + 0(6)] 



(76) 



As discussed in [12|, |38|, |46[ , this contribution to the wave- 
form will be very difficult to detect as it is essentially DC 
in character corresponding to a steadily-growing part of 
the waveform with no dependence on the orbital phase. It 
does, however, build up during the entire inspiral which 
leads to its magnitude being comparable with the OPN 
quadrupolar term despite arising formally as a 2.5PN 
contribution. See [13] for a discussion on strategies for 
detecting the memory piece of the waveform. 



B. Spin- weighted spherical harmonic components 



Examining Eq. (|T9)) . seems to show that the individ- 
ual spin-weighted spherical harmonic components could 
obtain contributions from both mass and current radia- 
tive multipoles. It turns out for the non-spinning case, 
however, that the mass (current) multipoles only con- 
tribute to components with i + m even (odd) . Because 
of this separation, and since C/^"* = 
and U^'"* = (-1)™F^.-'", it follows from Eq. ^ that 



(-l)™/i' 



{£ -\- m even) 



|(-l)™+i/i^"* (e + modd) 
which simply reduces to 



(77) 



(78) 



Substituting the source mass quadrupole (|47p into the 
expression for the radiative mass quadrupole (|23p . taking 
the appropriate time derivatives, substituting the equa- 
tions of motion (j40p . and evaluating the hereditary in- 
tegrals using the techniques described in Sec. IIV Al the 
dominant component of the waveform is given to 3PN 
order as 
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where the constant ro has been ehminating from the expression by using Eq. (|68p . The remaining components are 
given by 
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In Appendix [Cl spin-weighted harmonic components 
for even (odd) £ + m are given to IPN (OPN) order for 
arbitrary £. Again, the m < components are given 
by Cl]). 



C. Absorbing amplitude terms into the phase 

The ln(a;/a;o) terms that appear in the spin- weighted 
spherical harmonic components can be absorbed into 
a redefinition of the phase by introducing an auxihary 
phase variable ij) = (f> -\- 5. Since the In (x/xq) terms first 
enter at 1.5PN order, it is straightforward to show that 
choosing 0, El 



TO \Xq 



(117) 



where M = I (the mass monopole of the source), will 
eliminate the In (x/xq) terms from the components. This 
follows from 

= /i^™e-""'^[l - imS - ]W5^ + 0(a;9/2)], 

where h^™ is /i^™ omitting the ln(x/a::o) terms. Further- 
more, since the orbital phase as a function of frequency 
goes as 



1 



0(e). 



(118) 



at leading order, the In (x/xq) terms, which were 1.5PN, 
2.5PN, and 3PN order in the amplitude terms, now ap- 
pear as phase corrections at relative order 4PN, 5PN, and 
5.5PN. As these terms are beyond the order to which the 
orbital phase evolution is known (3.5PN order), it can be 
argued that these terms can be ignored. Note that the 
choices of xq in Eq. and S in Eq. (|117p are not unique 



as other amplitude terms can be absorbed into the phase 
(e.g. see [13); these choices were made to gather all log- 
arithmic terms into one term, as well as to simplify the 
waveform [isj . 

In order to recover the 2.5PN polarization waveforms 
in [l3, [11] from the components listed above, it is nec- 
essary to substitute ln(2:/a:;o) — *■ 0. After substituting 
the coefficients into Eq. (|lip and truncating the sum at 
2.5PN order, the result must be evaluated at Q = i, 
$ = 7r/2. Furthermore, there is an overall sign differ- 
ence due to a different choice of the polarization triad 
{N,P,Q). 

The polarization waveforms have been computed in the 
limit J/ — > in Q using black hole perturbation theory. 
In order to compare with their results, it is necessary to 
substitute 



Inl^ 

Xo 



17 2, „ 

h - n2, 

18 3 



into the /i^™ listed above as 0] makes a different choice in 
redefining the phase variable (and works in Schwarzschild 
coordinates as opposed to harmonic coordinates). After 
this substitution, and setting = and 5m /m — — 1, 
it is found that the /i^™ above agree with the results of 
0,[6i.2 



V. CONCLUSIONS 

It has been shown that the spin- weighted spherical har- 
monic components of the waveform can be computed to 
higher post-Newtonian order by computing them directly 
from the radiative multipolc moments rather than by pro- 
jecting them from the full polarization waveforms. In 



^ There is a sign difference in C,+ ^ between [6^ and 01 . The results 
presented here agree with the sign of [66| . 
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particular, this allows the dominant h?'^ component to be 
computed to 3PN order. Since numerical simulations can 
compute the waveform over the entire sphere, it is pos- 
sible to compare the spin-weighted spherical harmonic 
components of the waveform from the simulation with 
those predicted by a quasi-adiabatic post-Newtonian in- 
spiral. Thus, by examining the h?'^ component, it can 
be determined whether or not the 3PN contribution im- 
proves the agreement between the post-Newtonian wave- 
form and the numerical waveform. If significant improve- 
ment is found, it would suggest that it would be worth the 
effort of computing the full 3PN waveform in order to im- 
prove detection of marginal si gnal s 2^ 2^ 2?! a s well as 
improve parameter extraction [28L I29l.l30l.l3ll.l32l. Issl . Isl ] . 

In [23! a high-accuracy comparison is made between 
post-Newtonian generated waveforms and waveforms 
from a numerical simulation of 15 orbits of an inspiral 
of an equal-mass non-spinning binary black hole system. 
For this case, [23,] finds that the 3PN contributions to the 
amplitude of the /i^^ mode improve the accuracy with 
respect to the numerical waveforms. This suggests that 
for accurate parameter estimation, it may be desirable 
to compute the full 3PN amplitude for the polarization 
waveforms. 

For an equal-mass, non-spinning binary, only the 
source current octupole needs to be computed to have 
the full polarization waveform. For a non-spinning bi- 
nary with an arbitrary mass ratio, much more effort is 
required as explained at the end of [12| . But for com- 
parison with numerical simulations, post-Newtonian the- 
orists should keep in mind that extending the PN order 
of a given radiative multipole moment, will result in cor- 
responding improvements in the spin-weighted spherical 
harmonic components of the waveform, and these correc- 
tions will be of interest even if all the corrections needed 
to improve the polarization waveforms to the next order 
have yet to be computed. 
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APPENDIX A: IPN SOURCE MASS MOMENTS 



The source mass multipole moments for a system of 
N (nonrotating) compact point-masses to IPN order is 
given by m S ill 



N 



II = 



2(2£ + 3)c2 dt^ 



4(2^+1) 



(^+ l)(2^ + 3)c2 dt 



(Al) 



where m^i, j/a, and va are the mass, position, and veloc- 
ity respectively of the A-ih point mass, and 



i v\ G \ - TTlB 



f ^ ^ n n 



2 c2 c2 ^ \yA - VB 



.... ^^2) 
Restricting to the case of two bodies in a quasi-circular 

orbit, and transforming to the center-of-mass frame using 



y-i = !!^f[l-HO(e2)], 

TO 

2/2 = -^f[l + 0(62)] 
TO ■' 



(A3) 
(A4) 



I would like to thank Michael Boyle, Gregory Cook, 
Abdul Mroue and Saul Teukolsky for helpful discussions 



and using the relative equation of motion (|34p to elimi- 
nate the relative acceleration from the derivatives yields 



Xr, = vifi ■ 



hp- + 6£ 9 
2(£-hl)(2£-h3) 



X<L> 



£{£ + 9) , . .r 



2(£+l)(2^ + 3) 



r 



where 



where se = (m^-|-TO2)/TO^ and dg = {ml — 1112)/ , which 
themselves can be rewritten as polynomials in as 



/2 r 



{to, AM} 



{to, Sfe(z^)} for £ even 
{—6m, dkii^)} for £ odd 



(A6) 



Se 



fe=l 



£-k-l\ f£-k 
k-l 



k=0 



£-k-l 



(A7) 
(A8) 
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APPENDIX B: OPN SOURCE CURRENT 
MOMENTS 

The source current multipole moments for a system of 
N (nonspinning) compact point-masses is given to OPN 
order by [H, [ssj 



tributions from the radiative mass muhipoles. To IPN 
order, the radiative moments are given simply i time 
derivatives of the source moments. For arbitrary ^, the 
mass moment is given by Eq. (|A5p . In evaluating the 
spin-weighted spherical harmonic components from the 
source mass moments, the following identity is useful: 



N 



Jl 



"i^e Ha 



0{e) 



(Bl) 



Restricting to the case of two bodies in a quasi-circular 
orbit, and transforming to the center-of-mass frame as for 
the mass moments in Appendix [Aj it is straightforward 
to show that the source current multipole moments to 
leading order are given by 



Jl = —9e{^)L<i,XL-i> + 0(e), 
m 



(B2) 



r riLyL 



2'' 



(CI) 



where </> is the orbital phase (as opposed to an angle on 
the sphere). Since 4> = Lot, 



Y 



im 



=(-zmu;)'^-r^'"(|,</.). (C2) 



where 



Therefore, it can be shown that 



{m,5fc(zy)} = 



{-Sm,dk{iy)} 
{to, Sk{iy)} 



for 
for 



even 
odd 



(B3) 



eii-i) 



(f,</>), (C3) 



APPENDIX C: IPN SPIN- WEIGHTED 
SPHERICAL HARMONIC COMPONENTS 



using that f and lu are of 2.5 order. Finally, using that 



1. £ + m even 

For even £ + m (for non-spinning binaries), the spin- 
weighted spherical harmonic components are due to con- 



//2 



(C4) 



it can be shown that 



x£ Stt I {£+!){£ + 2) Gum tli^ri-mC^ \ \t i \ 



3 



-^fiM - ^ ^(/f{)(^/|3) ./£+i('^) + 0(63/2)| (for £ + TO even). 



(C5) 



where to and fk {i') are given by (jA6 



' -I- m odd 



current moment is given by Eq. ljB2[) . In evaluating the 
spin-weighted spherical harmonic components from the 
source current moments, the following identity is useful: 



For odd £ + m (for nonspinning binaries), the spin- 
weighted spherical harmonic components are due to con- 
tributions from the radiative current multipoles. To OPN 
order, the radiative moments are given simply £ time 
derivatives of the source moments. For arbitrary £, the 



6 ■ UL .yr - 1 / (2^+i)(^^-TO^) 

(C6) 

Following the same steps as above, it can be shown 
that 
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' = (2m)!! V m - m +m ~ d ^ b'^) (for^+-odd), 



where rh and (7^ (z^) are given by (|B3I 
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